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Abstract 


We consider a classical fermion and a classical scalar, propagating on two different kinds 
of 4-dimensional diffeomorphism breaking gravity backgrounds, and we derive the one-loop 
effective dispersion relation for matter, after integrating out gravitons. One gravity model 
involves quadratic divergences at one-loop, as in Einstein gravity, and the other model is the 
z = 3 non-projectable Horava-Lifshitz gravity, which involves logarithmic divergences only. 
Although these two models behave differently in the UV, the IR phenomenology for matter 
fields is comparable: (i) for generic values for the parameters, both models identify 10 10 GeV 
as the characteristic scale above which they are not consistent with current upper bounds on 
Lorentz symmetry violation; (ii) for both models, there is always a fine-tuning of parameters 
which allows the cancellation of the indicator for Lorentz symmetry violation. 


1 Introduction 

The high-energy behaviour of gravity is still not well understood, and several directions 
have been developed in order to build a satisfactory theoretical and phenomenological frame¬ 
work. One possibility is to consider gravity as an effective model jjLj and study the implications 
of its quantisation at energies well below the Planck mass. Another possibility is to modify 
the ultraviolet (UV) behaviour of gravity, as can consistently be done, for example, with 
Horava-Lifshitz gravity |2], We consider here two modifications of Einstein gravity which 
are invariant under the reduced symmetry of 3-dimensional diffeomorphisms, in the situation 
where a specihc coordinate frame is preferred, and study the consequences of this violation 
of local Lorentz symmetry on the low energy physics. In order to test the validity of these 
modified gravity models, we couple them to classical matter fields and study the generation, 
through quantum gravity corrections, of Lorentz-violating features in the matter dispersion 
relations. These two models are 

• A modification of Einstein gravity, which preserves only 3-dimensional space diffeo¬ 
morphisms, as in Horava-Lifshitz (HL) gravity, but which does not involve higher-order 
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space derivatives or anisotropic scaling between space and time, so that one-loop correc¬ 
tions are quadratically divergent. This model is consistent with the effective approach 
to gravity; 

• An extension of the first model, with an improved UV behaviour. This is achieved 
by the introduction of higher-order space derivatives, leading to less divergent loop 
integrals, and which can consistently be implemented if we allow for an anisotropic 
scaling between space and time. The extension we consider is the non-projectable 
version of HL gravity (npHL) |3] for z = 3, and involves logarithmic divergences only. 

In both cases, gravity is naturally described in terms of the Arnowitt-Deser-Misner (ADM) 
decomposition of the metric, which naturally exhibits a space-time foliation. Also, the gravi¬ 
ton has 3 degrees of freedom (dof), with a scalar dof appearing as a consequence of the 
breaking of 4-dimensional diffeomorphisms [4], 

We couple these models to classical matter and derive the one-loop effective dispersion 
relation seen by particles, after integrating out gravitons. One-loop quantum corrections to 
the dispersion relations of quantum matter fields coupled to HL gravity have been studied 
in [5], where the Authors derive the effective speed seen by a scalar held and an Abelian 
gauge held, and compare these to measure Lorentz-symmetry violation. In pjj, the non- 
projectable version of HL gravity is considered for the derivation of the effective matter 
Lagrangian. We are doing a similar study here, taking into account a classical scalar and 
fermionic backgrounds though, and we calculate the difference An 2 between the effective 
speeds of light seen by these two species. 

The hrst model involves quadratic divergences, so that we use a cut oh to calculate one- 
loop graviton loops, whereas the second model involves logarithmic divergences only, that 
we control with dimensional regularisation. We find that, although the two models behave 
differently in the UV, their phenomenological predictions are very similar in the sense that 
they predict the same order of magnitude for An 2 , if one keeps generic values for the different 
parameters. Therefore this phenomenology is not really improved by the introduction of 
higher-order space derivatives. Moreover, in both cases, we find that one can always impose 
An 2 = 0 if the parameters are fine-tuned. 

Studies of effective dispersion relations for Lifshitz-type models in flat space-time have 
shown the limitations of phenomenological viability of these models. This was hrst noted in 
[7], where the unnatural fine-tuning of bare parameters is shown in order to match the light 
cones seen by two different scalar particles interacting. Similar studies were done in [8], where 
the effective dispersion relation for interacting Lifshitz fermions is derived, in the case where 
flavour symmetry is broken, or in |9j, where the fermion effective dispersion relation is derived 
in a Lifshitz extension of Quantum Electrodynamics. Our aim here is to study similar features 
and, more precisely, how global Lorentz-symmetry is affected from a local symmetry breaking. 
A more complete study could involve the energy dependence of effective parameters obtained 
from a quantum modified gravity, in a Wilsonian renormalisation framework for example, 
as shown in ins for HL gravity. The latter approach is also used in m, where instead the 
scalar held is integrated out on classical metric background. 

Section 2 describes the models and discusses gauge freedom, as well as the actions for 
gravity and matter sectors. We describe the different simplifications which can be used for 
the one-loop integration of graviton, which is presented in section 3 and in the appendix. 
We integrate there successively the different components of graviton huctuations, for both 
models, where we keep on-shell the auxiliary helds which appear in the decomposition of the 
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graviton fluctuations. Finally, we discuss the phenomenology of the two models in section 4, 
based on the current upper bounds for Lorentz-symmetry violating parameters. 

2 Models 

In this article, Greek letters denote space-time indices and Latin letters space indices 
only. The signature of the metric is (—, +, +, +). 

2.1 The original actions 

We introduce here two modified gravity models, and their coupling to matter. 

• The first modified gravity action we consider is a modified Einstein-Hilbert action 
where new operators are allowed due to the reduced symmetry of the theory. Using the 
ADM decomposition and omitting the cosmological constant, it can be expressed as 

S G — Mp J dttd i xyfgN (AA, K lJ - XK 2 + R (3) + otarf ) , (1) 

where Mp = (IOttGat) - 1 and 

9 = | det(flfij-)| ( 2 ) 

Kij = ^(dtShi ~ D,N, - DjN t ) with /L,V ; 

K = K t ,g ij , R® = R ;jkl g ik </ , a t = d, In N . 

Gn and are respectively the Newton gravitational constant and the extrinsic curvature. 
The possibility to have A ^ 1 and a / 0 implies that 4-dimensional diffeomorphisms are 
broken to 3-dimensional diffeomorphisms, and the Hilbert-Einstein action can be recovered 
with A = 1 and a = 0. 

• The second model we consider is the non-projectable version of Horava-Lifshitz grav¬ 
ity m which allows the lapse function N to be a function of space and time. This implies that 
terms containing a* = di In N should be included in the action, in order to have a consistent 
dispersion relation for the scalar graviton [3], and avoids the well-known pathology present 
in the original HL gravity version |2]J. In this context, we impose space and time to scale 
anisotropically such that 

x —> bx while t —>• b z t , (3) 

and the choice z = 3 motivates us to introduce dimension 4 and dimension 6 operators, 
which become important in the ultraviolet (UV). As in [3j though, the time component is 
reparametrised so that the action is written in terms of the “physical” units. The npHL 
action can be written as [H] 

S HL = Ml j dttfxyfgN {KijK^ - XK 2 + + aa i( j (4) 

+F 1 R ij R ij + F 2 (A (3) ) 2 + FsR^Vitf + F^Aa* 

+Si(ViRjk) 2 + S 2 (\7iRW) 2 + 5 3 (Ai2®V i a i ) + S 4 (a,AV)} , 
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where Fi = ( fi/M‘j 1L ) and S) = (. Si/Mfj L ) with Mhl being the Horava-Lifshitz scale, and /,: 
and Sj are dimensionless coupling constants associated with operators of dimension 4 and 6 , 
respectively. 

• Finally, for the matter sector, we will consider complex scalar and fermion fields mini¬ 
mally coupled to the gravity models ([I]) and fl3J). The complex scalar field action is 

S s = ~ J dtd 3 x yfgNg^d^dvff)* , (5) 

and the fermion action is 

S,=~ j dtd?x j [h a e^V^ - e'JVrfb"V>] . ( 6 ) 

where 

e = det(e/) = y/gN (7) 

+ and Vy 0 = 

= ^w^ a pa a P and a ap = ^[ 7 “, 7 ^] , 

— e a {d)fiG\/3 ~ ^ Xypcrp) — e a(-^V e A/?) j 

^a /3 being the spin connection. 

2.2 Gauge invariance and degrees of freedom 

The 4-dimensional diffeomorphisms are explicitly broken for both gravity models for A 7 ^ 
1, a ^ 0, Fj ^ 0 and Si 7 ^ 0. Instead, these models are invariant under 3-dimensional 
diffeomorphisms 

st = m ( 8 ) 

5x l = £ l (£, x) 

ddij d)i£,j T dj^i T £ dkQij T f ijij 

SNi = d i £ k N k + Z k d k N i + &g ij + fN i + fN i 

5N = £ k d k N + fN + fN . 

Because of 4-dimensional diffeomorphism breaking, a third physical dof is present in both 
models, and a way of counting the number of dof of the theory is given in [12], where the 
Authors use a Hamiltonian description of gauge theories. They show that the number of 
primary constraints to take into account is the number of gauge functions plus the number of 
their time derivatives, in the situation where these gauge functions depend on space and time. 
This is because gauge functions and their time derivatives must be considered independently, 
when defining a boundary condition for the evolution of gauge fields. In our case, we have 10 
independent metric components (N, Ni, g t j) and we see from the gauge transformations (151) 
that the functions count twice since they appear with their time derivative, while / counts 
once only because it depends on t only. The total number of dof is therefore 10 —(2x3+1) = 3. 
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In order to study the one-loop quantum corrections to the matter sectors, we expand the 
metric g, w and, consequently, around a flat background: 

9[iv 9/iu T hfiD (9) 

g lw = rT -hT + hP x h v x + -- - , 

< = ^ + 1 -h“- 1 -h^ + --- 
e\ = K ~ \K + + ■■■ , 

where dots represent higher orders in fluctuations. Using the following relations 

9^ = efefVap , ( 10 ) 

p a P V — n a p a p fl — a 11 

e yS ~ 'Ip i e a — 9u i 

g^ v = —N 2 d 2 t + gij(dx l + N l dt)(dx J + N^dt) , 

we can then express the different ADM components in terms of their fluctuations 


N 

= 1 + 71 

(11) 

Ni 

= rii 


9ij 

A hjj . 



The fluctuations n l and h t] can additionally be decomposed into their different spin compo¬ 
nents as: 

rii = nf + dip , (12) 

d 2 ^j B+ S -fh, (13) 

where H^ is a transverse-traceless tensor, nf and W t are transverse vectors and 5, h and 
p are scalar fields, h being the trace of hij. Nevertheless, making use of the gauge freedom 
shown in eq. (j8j), a natural gauge choice is to set W % and scalar B to zero. Consequently, 
eq. (H3]) becomes 

r 

ha = + -fh , (14) 

where and h are the 3 physical degrees of freedom present in the theory, while n, nf and 
p are auxiliary fields. 


hij = + (diWj + djWi) + did. 


jij 


2.3 Expanding the actions 


Since we are interested in one-loop corrections, it is enough to expand the actions up to 
quadratic order in the ADM field fluctuations. The flat space metric is Sij, such that all the 
spatial indices are lowered, i.e. h—> h %3 . We have then 


V~9 

r fc 


1 + —h + — (/r — 2 hijhij) + • • • 
z o 


a. 


2 


d,h k j I hudihj k T hijdihki ,, hjh,i 


(15) 


where h = ha, and dots represent higher orders in fluctuations which can be omitted for the 
one-loop calculation. 
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2.3.1 Matter sector 


We explain here the construction of the relevant matter actions for scalars and fermions, 
and then describe the ansatz taken for these external fields. 

• Expanding the scalar action (jSJ) up to quadratic order in the graviton field fluctuations, 
the scalar action becomes 


= 


dtd 3 x 


h hn l.o ,o 
l + n+ - + — + -(h - 2 h. 


hm 


ij) 


- 00 * + d k <pd k (f)* 


(16) 


2 n 00 * + 2rij0<9j0* — hijd^dj^ + 2 ( —-nn* — rijhij 


+ (nh — u)(j>q !>* + ( huhij — riiUj — nh. 


hhj 




di4>dj4> 


3 7 r ) 


where the first line can only generate Lorentz-symmetric contributions and therefore will be 
omitted. Moreover, as pointed out in [5j, quadratic terms in the metric fluctuations will 
only contribute to our results when the graviton fields are contracted among themselves. 
Therefore, for any tensor 70- quadratic in the graviton field, because of rotational invariance 
in space, one can use the following simplification: Tijdi^djcj)* —$■ (Ta/ 3)9^09^0*. In addition, 
linear terms in the metric perturbations can also be omitted, as they lead to quartic matter 
self interactions after completing the square (we are interested in corrections to the kinetic 
terms only). Finally, terms of the form hn ,, nn t or rijhij do not contribute since they have 
to be contracted with another vector metric fluctuations, leading to terms which are at least 
cubic in fluctuations. The relevant part of the action, containing only terms which generate 
Lorentz-violating contributions, is then 

<S^ 2 ) = — J dtd 3 x 


(nh — n 2 ) 00 * 


J (y - n\ - nh - y ] 


(17) 


To simplify even more the expression above, we can write 

(nh — tt, 2 )00* = — (nh — n 2 )^0(9^0* + (nh — n 2 )<9fc0<9fc0* , (18) 


and since the first term on the right-hand side of the expression above is Lorentz-symmetric, 
we only need to take into account the second one, then (fI7|) becomes: 



dtd 3 x 


Kj ~ n i —2>n 


+ 2 nh — 


h 2 


<9fc0<9 fc 0* 


(19) 


For the fermion sector, we find 


o(2) _ _ l 

b f ~ 9 


/■ f 

h hn 1 9 0s 

/ dtd 3 x < 

I + 77 + - + — + -(h - 2 hij) 


07^0 


+ ( nS ij~2 hij ) (^c^0 + 0{7n r j i) }0) +0{7 / ",rd ) }0 


1 

+ 2 U 


0 (ri'S - 7°W + { 7 u rj^} - {7°, rf } }) 0 + 0{y M , rj t 2) }0 
hn l 3 

T 




77.77,1 


'-rijhij - ^ 1 0 ( 77 ^ - 7° c?) 0 + ^ 0 { 7 /i , rj t 1 ) } 0 + 


1 9 T ntd-, 03, , hha n(hSa — ha) 

+ -n^a^ + -huh,, - -y + — 'ii 


\n { nj ) 07^0 


( 20 ) 
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As in the scalar case, the terms on the first line would only contribute to Lorentz-symmetric 
corrections, hence such terms will be omitted. For the same reasons explained above, linear 
terms in the graviton fields as well as terms proportional to 07 \ do 0 or 0 70 <9* 0 do not con¬ 
tribute to our calculation and therefore will also be omitted from now on. Terms containing 
the anticommutator of a 7 matrix and the spin connection will not contribute as well, since 
they involve the contraction of a symmetric tensor with an antisymmetric one. Finally, as 
discussed above for the scalar action, when T V} is a given quadratic term in graviton, we make 

_ ^_ y _ ^_ y 

the replacement 70 - 07 j< 9 j 0 —>■ (Tk/3)(07fc dk 0)- Thus, the Lorentz-violating fermion action 
we consider becomes 

Sf ] = - l - f dtd 3 x j^n-07°o?0 + ^ Qh 2 - - +nh - n 07^0 j , (21) 

which, after removing Lorentz-symmetric terms, reduces to 


o( 2) 

W 


/ dtd 3 x 

-hi 

h 2 

1 2 2 ' 

- n ,■ H —nh 

} 

4 l o 

6 

2 * 3 


in k dk 0 


( 22 ) 


• After reducing the matter action to their simplest form (TT9l) and (l 22 ]h we then consider 
a plane-wave as an ansatz for the external fields 


00*0 = 0 o exp (-iffx,,) 

0(0) = 0o exp {-iq^x^) , (23) 


so that the quantities C00C00* and 0 i 7 fcC 00 become constants and, therefore, will be re¬ 
spectively replaced by and (— 2 f)[ 0 o (7 • <y)0o]. With these classical background matter 

held configurations and using the held decompositions (1T2|) and <!» the scalar ( 11911 and 
fermion (|22|) actions are respectively 


Sf> = 

1 

[ dtd 3 x 

Hi - 

h 2 

3 

l 3 

6 

o( 2 ) 

b f - 

1 

[ dtd 3 x 

-Hi 

h 2 

2 J 

4 v 

12 


(nj ) 2 + pd 2 p — 3 n 2 + 2 nh 
2 


(nf) 2 1 2 

— +- 2 p d p+r h 


tf<&) , (24) 

[0o (7 • q)M ■ 


2.3.2 Gravity actions 

For the gravity actions, we expand (P and (QJ up to quadratic order in the metric fluc¬ 
tuations and make use of the metric decompositions (fl 2 |l and (fT4l) to obtain 


S { a ] = Ml 


dtd 3 x 
(3A - 1) 


~ d 2 )^ - ^njd 2 nj - (A - l)p{d 2 ) 2 p 


(25) 


12 


hd+h —— hd 2 h — and 2 n + —^-—p<9 2 h — -nd 2 h 

t 18 3 r 3 
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and 


c( 2 ) 

d hl 


(26) 


= Mp I dtd 3 x <J \h i3 [-d 2 + d 2 + F\(<9 2 ) 2 - S^d 2 ) 3 ] H %3 - ^njd 2 nj 


18 


3(3A — 1) 


d 2 -3 2 + (3 F 1 + 8F 2 )(«9 2 ) 2 - (3S 1 + 8S 2 )(d 2 ) 3 


n[ad 2 + F 4 (<9 2 ) 2 + S 4 (<9 2 ) 3 ]n - (A - l)p(«9 2 ) 2 p + (3A ^ ^ 2 j 


-pd 2 h 


- h[d 2 + F 3 (d 2 ) 2 + S 3 (d 2 ) 3 ] 


n 


Finally, because ghosts do not couple to the matter sector at tree level, one needs to consider 
at least two-loop corrections in order to have a non-vanishing contribution coming from 
interactions with ghost fields. Therefore, in the present work, ghosts can be omitted. 


3 One-loop matter effective actions 

We integrate here the metric fluctuations in order to obtain the effective matter kinetic 
terms. For both gravity models, we impose the Hamiltonian and momentum constraints in 
the path integral, which consists in keeping auxiliary fields on-shell. This approach is used 
in [ 6 ] and implies that no conformal instability arises in our calculations |l3j . It is known 
in perturbative quantum gravity that, when introducing an irreducible decomposition for 
the metric, some of the components have a propagator with the wrong sign, leading to a 
potential problem in defining the partition function. This unstable mode can be traced down 
to a conformal factor, but with our gauge choice it would arise from the integration of the 
auxiliary fields in the metric decomposition, which does not happen here. 

This problem of conformal instability can be understood as an artifact arising from per¬ 
turbative expansion [14], but can also be avoided by making an analytical continuation of the 
metric components to imaginary values, simultaneously with the Wick rotation [15]. We note 
here another non-trivial connection between the Wick rotation and quantum gravity, in the 
context of time-dependent bosonic strings. As shown in [TB] for a specific string configura¬ 
tion, which satisfies conformal invariance non-perturbatively in a', a Wick rotation in target 
space implies a phase factor for the overall string partition function. This phase becomes real 
for specific space-time dimensions only, and 4 appears to be the lowest dimension for which 
the string partition function remains real, after the Wick rotation. 

3.1 Model I: modified Einstein-Hilbert gravity 

We study here the model ([I]), for which one-loop quantum corrections are quadratically 
divergent, and will therefore be regularised with a cut off. After expansion of the matter and 
gravity sectors in terms of the metric fluctuations, the actions that we are interested in here 
are (l24|) and fl25jh 

3.1.1 Constraints 

Fluctuations in the shift vector (nf,p) and the lapse function (n) are auxiliary fields 
and are therefore not propagating, such that one can use their equations of motion (the 
momentum and Hamiltonian constraints) that we substitute back into the action. 
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Varying the actions with respect to n gives us the following constraint 


[-2 aM P d + 2{p (j> 0 )]n = - 


M P d + (p 0 0 ) + -[^ 0(7 • q)M 


h 


(27) 


whereas variations with respect to the scalar and transverse parts of the shift vector lead to 


2Mp(l - \)d 2 ~-[(p (j), 0 ) + -[^ 0(7 • q)i> 0 ] 


d 2 p = 


M 2 P { 3A- 1) 


- M 2 P d 2 + 


iP 0o) + 7 [^ 0(7 -fill’d. 


ni = 0 


d 2 h , (28) 


(29) 


When the last constraint is put back into the actions, all contributions coming from nf 
disappear and from now on such actions will only depend on the tensor and scalar components 
of the metric. On the other hand, since the auxiliary scalar fields n and p appear mixed with 
the scalar graviton h, before substituting the constraints ((27 J and (128]) back into the actions, 
we can expand them in terms of the matter contributions to find 


n = — - 


P = 


3a 1 + 
(3A - 1) 
6 (A - 1) [ 


(pVo) f a + 1 


Ml 


a 


(d z r + 


[^0(7 

2 Ml 


{& 


, 23—1 


h + ■ ■ ■ , 




(30) 

(31) 


where dots represent higher-order terms in the matter fields. The other equations of motion 
fed back into the actions lead to 


s[ 2) = 


dtd 3 x <! -Hi. 




2 (<9 2 - d 2 ) - |(p®^) - ^ 0 ( 7 '9)^o] 


H 


-h 


M 2 


-Xd 2 + 


a 


a 


d 2 + 


fa 2 + 4a + 2 X 2 d 2 


+ 


cr 


6 d 2 


+ 


[^ 0 ( 7 ' 0 )^ 0 } 13a + 8 X 2 d 2 t 


4ck 


+ 8 d 2 


(32) 


where 


3A - 1 


X = — . (33) 

For a consistent propagation of the scalar graviton h, one needs X > 0, such that the allowed 
values for A are 

A < 1/3 or A > 1 . (34) 

Finally, from the action (1321) we find the following dispersion relation for h 


lu 2 = 


A — 1 
3A - 1 


2 — a 


a 


k 2 , 


(35) 


which shows a consistent propagation for 0 < a < 2 for the allowed values of A given in 

eq.flMD. 
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3.1.2 Loop integration 


We give in the Appendix the details of the integration over graviton. 


Spin-2 component 

The integration over the spin-2 component gives 


exp 


M 2 P 


^(pVo) + -q)M 


m 

2(2vr) : 


■A 2 H- 


where 5(0) is the space-time volume and dots represent either field-independent terms or 
higher orders in (p 2 ^) and [0o(7 '<p)0o]- 


Spin-0 component 

The integration over the spin-0 component gives 


exp 


M 2 P L 


a 2 T 4 a f 2, n o N 3a T 8 


2a 2 


m 0 ) + 


8 a 


[-00 (7 • q)M 


( (0>o) , [0o(7-g)0o] 
4 V 3 4 


r ^(0)A 2 . 

2(2t r) 2 


where 


Y = 


a(A — 1) 


(2 — a)(3A — 1) 


(36) 


(37) 


3.1.3 Total Lorentz-violating contributions 

Considering the results obtained in the previous sections, we add here all relevant contri¬ 
butions and present the total Lorentz-violating corrections for both scalar and fermion fields. 
We also note that 


(pWo) = 


5(0) 


dtcl 3 x d k (j)d k (f)* and [0o(7 ' <h)0o] = 


5(0) 


dtd 3 x i\)i^ k d k i\) . (38) 


Considering the results (1671) and CD, the total contributions to the matter fields are 
given by the following Lagrangians 


A 2 


2(2vr) 2 M 2 


in the scalar case, and 


p L 


A 2 


u 


(X 2 

a 2 + 4a + 2 A 

[s 

2 

{ 6 

a 2 ) 


d k (j)d 


(39) 


T 


( 3a + 8 

x 2 y 

_2 

8 

^ a 

2 )\ 


027aA0 , 


(40) 


in the fermion case, where X and Y are defined in eq. (]33l) and eq. f|73]l respectively. 
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3.2 Model II: non-projectable HL gravity 

We turn here to the non-projectable version of Horava-Lifshitz gravity (J4J) , for which the 
relevant actions are given by eqs. (l24p and (126)) . In the absence of matter, one-loop quantum 
corrections for this model are logarithmically divergent, but in the presence of dynamical 
matter, quadratic divergences arise from the coupling gravity-bosonic matter M. In our 
case, matter cannot induce further divergences compared to the single gravity case, since 
it is classical and thus does not involve any new loop momentum. The use of Hamiltonian 
and momentum constraints generate an artificial quartic divergence though, due to the in¬ 
troduction of additional space derivatives in the decompositions (jl2M13j) of the graviton: as 
can be seen from the action (|44l) below, the presence of matter is then accompanied with the 
derivative operator <9 2 /(<9) 2 , arising from the coupling between p and h in the gravity sector. 
This divergence is thus a gauge artifact, on which we will come back in section 4, and that 
we disregard in our calculations. 

As a consequence, we use dimensional regularisation, where the naively power-law diver¬ 
gent integrals actually vanish Ed- We note that the vanishing or finiteness of a regularised 
integral which otherwise would naively be divergent is explained pedagogically in [TS]: in the 
regularised integral, divergences associated to different regions of the domain of integration 
cancel each other, such that the integral is finite when the regulator is removed. 


3.2.1 Constraints 

The constraints (128]) and (1291) obtained from the variation of the action with respect to 
the shift vector components (nf,p) are the same as in the previous modified gravity model, 
since the dimension 4 and dimension 6 operators which are added do not depend on these 
components. The additional contributions to the lapse function fluctuations n lead to the 
new constraint 


[—2MpT>2 + 2(p 2 0 2 )]n = - 


M P V i + (p </> 0 ) + - fata (7' 0)M 


h, 


which can be written as 
1 


n = — - 


Pi (p 2 ^) 1 

V 2 M 2 P V 2 


1 + 


Pi 

Vo 


+ 


# 0(7 • q)iPo] 


2 Afl lTUW ,/TUJ D, 


/& + ■■■ , 


where dots represent higher-order terms in [^ 0(7 • Qj^o] and (fPcj) o)> and 

Pr = [ d 2 + F 3 (D 2 ) 2 + S'sPa 2 ) 3 ] , 

V 2 = [ad 2 + F 4 (<9 2 ) 2 + s 4 {d 2 ) 3 } . 


(41) 


(42) 


(43) 


Using the constraints (|28]h (l29|l and (1421) to rewrite the original actions (124)) and fl26|) . we 
arrive at the following action, which only depends on the physical metric fluctuations H, t j 
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and h, 

o( 2 ) _ 

— 


dtd 3 x < -Hi, 


■-IJ 


M(_ 9? + a 2 + n(a 2 ) 2 - ^o 2 ) 3 ) -1(^0?) - \[Mi■ «)*] 


-/l 

2 


^ f~Xd 2 - a 2 + (3F, + 8F 2 )(<9 2 ) 2 - (35, + 85 2 )(<9 2 ) 3 + 2 f 2 


1+4, ^J + H^J + irl 


(pV|) I, . . m\ . 001 

9 


[^o(7 ■ g^o] A3 A Pi\ X 2 ^ 

9 U + lx>J 8 a 2 


h 


Hi:i 


(44) 


3.2.2 Loop integration 

We give in the Appendix the details of the integration over graviton, which is done using 
dimensional regularisation, with d = 3 — e. 


Spin-2 component 

The integration over the spin-2 component gives 


exp 


Ml L 


^(pVo) + ^ 0(7 ■ q)i>o] 


m it 

(2vr)2^5^ e 


+ ... 


(45) 


Spin-0 component 

The integration over the spin-0 component gives 


exp 


Ml 


2 \ c, 


(pVo) + 7 ^ 0 ( 7 ' 9 )^ 0 ]) + 

(0) \ 4 J J c (2) 


c, 


( 1 ) 


(2(pV3) + [^ 0(7 • q)M) 


5(0) // 

(2vr) 2 v / X e 


+ 


where Cg , with n — 0,1, 2, are given in the Appendix. 


(46) 


3.2.3 Total Lorentz-violating contributions 

In the Appendix appears the following integral 

1(A) = f d “ k 1 


( 27r ) d kWk 2 + A ’ 
which, with dimensional regularisation {d = 3 — e), gives 

I(A) = 2^7 + ' 


(47) 


(48) 
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In the limit e —» 0, we obtain then 


4 i(a) =i ’ (49) 

such that 

X(A) = 2^ ln (^) ’ (50) 

where /x 0 is a mass scale. To choose /x 0 accordingly, we calculate 1(A) using a cut off A in 
dimension d = 3 and find: 


J(A) 


1 , / A + VA 2 + A 


2 vr 2 


In 


y/A 


( 51 ) 


From the form of A in (1751) and (1881) . we note that A can be written as dMl Ll where d repre¬ 
sents a dimensionless constant of order 1 for each of the different cases. Then, expanding (l5Tj) 
for A M hl , we find 

x w = & M (£r) ■ (52) 

where finite terms were omitted in the expression above. Comparing eq. (l50l) with the result 
above, we naturally choose /xo = Mhl and n — A. 

Considering now the results obtained above, we write the total Lorentz-violating contri¬ 
butions for both scalar and fermion fields. Using the relations (138|) and assuming (1501) with 
/io = Mhl , the total corrections for scalar and fermion fields are, respectively, 


1 

2 ( 2^)2 

1 

2 ( 2^)2 


4 1 11 

3 vW ^ \J ATg 0) 




Mhl ]n (M% l \ 
M 2 P V A 2 ) 


d k (f)d k (j)* , 


11 3 1 

^VW\ 8 ^/x4 0) 



Mhl , 

m 2 v a 2 ; 


'i/ji'fkdk'ip . 


(53) 


3.3 Non-minimal coupling 

The models studied here involve minimally coupled matter, and one could ask what the 
effect of a non-minimal coupling would have on the effective dispersion relations. We show 
here that these would not affect our results. 

Given that we impose all the terms in the action to be at most of mass dimension 6, and 
that the scalar held is dimensionless for z — 3 in d = 3 space dimensions, its non-minimal 
coupling to gravity would be for example of the form 


M (3) + 6 ( R m ) 2 + 6(^)1 (#*) . 


(54) 


where G have the appropriate mass dimension for the terms inside the square bracket to be 
of dimension 6. Similarly, for the fermion of mass dimension 3/2, we could have 

[CiR (3) + C 2 (diCi 1 )] (-0-0) • (55) 
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In both cases, the space derivatives appearing in R (:i> and a* generate space derivatives for 
matter fields, after integration by parts, and thus could naively contribute to the effective 
dispersion relation. But the latter is obtained with matter plane wave configurations, such 
that and -0-0 are constants, and these configurations would therefore not give additional 
contributions to the dispersion relation. A more general held configuration could be chosen 
of course, but which would also contribute to all the other terms calculated in this article, 
in such a way that the final dispersion relation would not change: the functional for matter 
holds obtained after integrating gravitons is unique, and the corresponding dispersion relation 
is obtained by plugging a plane wave solution. 

This conclusion is valid at one-loop though: non-minimal coupling would radiatively gen¬ 
erate terms which would modify the matter kinetic terms, with an impact on the dispersion 
relation at two loops and higher orders. 


3.4 Comments on regularisation 

All the integrals in this article have been calculated by first integrating over frequency 
and then over momentum. In most of the integrals, the integration over frequency is finite, 
and has been performed without any regularisation. We are then left with a 3-dimensional 
integration over momentum, which is regularised either with a cut oh or with dimensional 
regularisation. But there is also the situation where the integration over frequency is diver¬ 
gent, and we discuss here few details for both models. 


• Model I 

As discussed previously, this model involves quadratic divergences and therefore cannot be 
treated with dimensional regularisation, which sees only logarithmic divergences. A typical 
example where the integration over frequencies is finite is 


r d A k 1 
J (2vr) 4 k >2 _|_ z 2 k 2 


1 

47T 3 




du- 


k 2 


u 2 + z 2 k 2 


A 2 

8n 2 z 


(56) 


where z is a positive dimensionless constant. 

But we also have integrals for which the integration over frequencies is divergent, and these 
are calculated with the same cut oh as the one used for momentum integration in the other 
integrals. The integral is of the form 


d 4 k 


u 


(2vr) 4 k 2 {u 2 + Z 2p) 


4-tt 3 

1 


2 7,2 


4-tt 3 

A 2 

8n 2 z 


,°° ,+A 2 

( dk / du -— 

0 J A U 2 + z 2 k 2 

pOO /»+A 

f dk du 1 -— 

o J- a u 2 + z 2 k 2 

poo 

I dk 


z 2 k 


A — zk arctan f — ) 
\zk 


(57) 


We note that the same result is obtained if one performs hrst the finite integration over 
momentum, and then uses the cut oh for frequency, whereas the commutativity of the order 
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of integration is not obvious when the integrals are divergent. 


• Model II 

The second model contains only logarithmic divergences, but artificial quartic divergences 
appear as a result of the graviton decomposition (fT3l) in terms of auxiliary fields. Since 
[a;] =3, we regularise the integral over frequency with the cut off A 3 


d A k 


LU 


(2 tt) 4 + z e k e ) 


-i poo p+A 3 

-—o / dk du 

4tt 3 ./o J_ A 3 


1 - 


3 fc 6 


1 


2vr 3 J 0 

A 4 


dk 


A 3 — z 3 k 3 arctan 


to 2 + z e k 6 
A 3 


(58) 


z 3 k 3 


8tt 2 z 


As explained at the beginning of section 3.2, this divergence is artificial and is thus omitted. 
The other integrals for this model involve a finite integration over frequency, and a logarithmi¬ 
cally divergent momentum integral. The latter is regularised with dimensional regularisation, 
but the same coefficient of the logarithmic divergence would be found with a cut off, as ex¬ 
plained in section 3.2.3. Therefore the Lorentz-symmetry violating terms generated in this 
model are unique (in the present gauge) and independent of the regularisation. 


4 Analysis 

The Lorentz-symmetry violating effect is measured by the product v 2 of group and phase 
velocities, which should be equal to 1 in the Lorentz-symmetric case. As noted in [5], the 
measurable deviation from the Lorentz-symmetric case is the difference |u 2 —v 2 |, which should 
be typically smaller than ICO 20 , according to the current upper bounds on Lorentz-symmetry 
violation pg. We note that |5] consider dynamical matter fields, instead of classical ones, 
which allows the cancellation of the above mentioned quartic divergence, as expected for a 
gauge artifact. Indeed, the graviton loop giving rise to this divergence is canceled by the 
equivalent matter loop, the difference of signs coming from the fact that the gravity spin- 
0 held component leading to the quartic divergence is auxiliary. In our case though, this 
cancellation does not occur because the matter loop is not present. This shows the non¬ 
trivial fact that a classical matter background can be consistently taken into account only 
if one removes this specific gauge artifact by hand. Other gauge artifacts are removed by 
taking the difference |u 2 — v 2 \, as we do below. 

The calculation of the traces obtained in the previous sections leads to the following 
kinetic terms 


+ v < f l \i>lkdkip (59) 

-dtMF+ivWydrfdrf* , 

where v™ — 1 + Sv™ and (w ™) 2 = 1 + ( 8 v ™) 2 and m — I, II for the first model (JTJ) and the 
second model (J4j) respectively. 
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4.1 Model I 


For the first model, ( Sv f) 2 and Svj are given by eqs. (l39j) and (HD)) , respectively, i.e. 


(Sv s ) ~ 


,I \2 


2(2vr) 2 


4 Y f X 2 ct 2 + 4a + 2 
3 + 2" 

1 


[l r 

/3a + 8 

x 2 y 

4 

^ a 

2 )\ 


A 2 

M 2 ’ 
A 2 


M 2 


( 60 ) 

(61) 


Subtracting the fermion contribution from the scalar one, we get the measurable departure 
from Lorentz symmetry 


IOf 2 - Of 2 | = |^(A,a)| , where 


F{ A, «) — — g + ^ 


(62) 


T— I 

1 

-< 

r (3A—i) 2 

a 2 — 4 

/ (2 — a)(3A — 1) 

"p 

1 

to 

a 2 


If one does not impose any specihc values for A and a, F( A, a) is generically of order 1, such 
that the current experimental bounds on Lorentz violation are satisfied with the difference 
(1621) if the graviton loop momentum is cut off by A < 10 10 GeV. 

Nevertheless, one can always find specihc values for A and a such that the difference (16211 
vanishes. An example is 


F(Ao,ao) = 0 for Ao — 0.332 and ao — 1.995 , (63) 

which are both in the range of allowed values for these parameters. This set of values has 
been chosen here because they are close to the boundary values A b = 1/3 and ot\, = 2, 
which suggest that quantum corrections point towards this specihc point in parameter space. 
It is interesting to note that \ = 1/3 is found as an IR hxed point for the Wilsonian 
renormalisation hows studied in and it would be interesting to relate this result to ours. 

4.2 Model II 


For the second model, from the results in ()53ll we can write (Sv^. 1 ) 2 and Svj 1 as 



Consequently, the physical quantity in which we are interested, i.e. the difference between 
these two contributions, is 


I Of ) 2 




Mhl ln (M% l A 
M 2 V A 2 ) 


(65) 
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If one assumes that A ~ Mp and Si, are of order 1, the bounds on Lorentz violation [20| 
are then satisfied for M HL < 10 10 GeV. Note that, by construction of HL gravity, M HL should 
also be large enough for higher-order space derivatives to be small compared to relativistic 
terms in the IR. The fact that quantum corrections imply an upper bound on M HIj is because 
we impose the UV regime, above Mhl , to dominate the loop integrals “early” enough, for 
quantum corrections not to be too large. 

However, similarly to what we find with the first model, the Lorentz-violating correction 
in (jggp can also be canceled if the coupling constants are chosen accordingly. For instance, 
one of the many ways in which this cancellation takes place is achieved by setting .Sj, c ^ to 
1, and taking A ~ 1.969, which is in the allowed regime for this parameter. 

4.3 Conclusion 

Our results lead to the following main points 

• If one wishes to conclude with generic values for the different parameters, then both 
models lead to the same order of magnitude 10 10 GeV for the typical scale above which 
the predicted Lorentz symmetry violation is too large. Therefore, although the second 
model has an improved behaviour in terms of UV divergences, the IR phenomenology 
is not really improved compared to the model without higher-order space derivatives; 

• If one accepts to fine-tune the different parameters, then there is always a range of 
values for these parameters, such that the effective maximum speed seen by particles is 
consistent with Special Relativity. Once again, this is valid for both models, and this 
feature is not a consequence of introducing higher-order space derivatives. 

Finally, the typical scale 10 10 GeV is consistent with other modified gravity model, as 
in [5]. The same characteristic scale is also found in [2Tj, where non-relativistic corrections 
to matter kinetic terms are calculated in the framework of the Covariant HL gravity [ 22] . 
We note that this scale also corresponds to the Higgs potential instability [23], which could 
be avoided by taking into account curvature effects in the calculation of the Higgs potential 
[24] . It would therefore be interesting to look for a stabilising mechanism in the framework 
of non-relativistic gravity models. 

A next step consists in deriving higher-order space derivatives for the fermion kinetic 
term, from a 4-dimensional diffeomorphism breaking gravity model, in order to generate 
fermion mass and flavour oscillations dynamically, as was shown in [25]. 


Appendix: loop calculations 

Model I 

Spin-2 component 

Gathering all the terms depending on LR,- in the actions 


and (1251) . we obtain 


Sf\H l3 ) = I dtd A x -H t] 


1 


/If 2 

^(«9 2 - d 2 


jj(pVo) - -?M>] 


H, 


ij , 


( 66 ) 
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which can be rewritten in terms of the Fourier transform Hij(k) of H t] (x) after a Wick 
rotation (t —)■ it, u —)■ —ico) as 


VHij exp <! iS[ 2 \H i:j ) \ ->• / VIL ri exp <! - 


1 - 

“W 2 H ” ih) 


A 


(i) 

Hu 


Hij(k 


V (> 


where .4.^ = 


Hi. 


+ |(pVg) + ^#0(7 • 9)^0] 


5{k\ + k 2 ) 


with k\ = cc 2 + k'i in Euclidean space. Then, considering the two components of Hij, we can 
perform the functional integration to find 


det 


M 2 2 1 

(2tt) 8 V "W' + + ) *(*i + fc 2> 


(67) 


= exp <j — Tr In 
1 


L(2 + 


5(^1 + fc 2 ) 


= exp 


M 2 L 


j(pVo) + tjW’oCt 


IVI dh±M I + 


where dots represent field-independent terms or higher orders in {p 2 ^ jj) and [if; 0 ( 7 -< 7 )^ 0 ] • The 
trace is calculated using a momentum cut off A and leads to 


Tr 


S(ki + k 2 ) 


d 4 M 4 fe fS{k 1 + k 2 )\ , 5(0) 2 


kf J J (27r) 

where 5(0) is the space-time volume. 


k 2 


5(ki + k 2 ) = 


2(2tt) : 


( 68 ) 


Spin-0 component 

Considering the terms which depend on the Fourier transform h in the action 
performing a Wick rotation, we find 


Vhex p iS 1 } 2 ' 


-> / Vh exp < - 


d A k\ d A k 2 1 

(2t r) 4 (2 tt) 4 2 


h{k 2 


A 


(i) 


h(ki 


where 


-4° = ^ mi (xu* + ) + (p ^) 


+ [^o(7-9)^o] 


CD 

X 2 uf 3a + 
8 kj 


X 2 bj\ a 2 + 4a + 2 

"6“ I 2 


a^ 


4a 


5(ki + k 2 ) . 


The evaluation of the functional integral leads to 


exp 


M 2 


“ 2 + 4 “ + V<® + ^[*(7 ■ 9>„] ) Tr 


2a 2 

X 2 ( (p 2 ^) [^ 0(7 • 


6 (fa + k 2 ) 


and 


(69) 


(70) 


+ 


Tr 


8« L, ” v ' "’“7 \Xo7 2 + a- 1 (2-a)Ai 2 

5(fci + /c 2 )<^ 2 


k'i (Xcuj + a _1 (2 — a)kj 


(71) 
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Both traces above are evaluated with the common cut off A for frequencies and wave vectors, 
and lead to the same result 


Tr 


d{k\ + k 2 )ujf 

k\ ^Xcj 2 + a~ 1 (2 -a)kfj 


= Tr 


S(ki + k 2 ) 

(Xuij + a-H2-a)k'i) 


v mT 

2 (2?r) 2 


,(72) 


where 


Y = 


a(A — 1) 


(2 — a)(3A — 1) 


(73) 


Model II 

Spin-2 component 

Comparing the terms which depend on Hij in the expressions (132|) and (l44]h it is clear that 
the only difference is on the appearance of higher-order space derivatives in the propagator 
of the tensor field in the npHL case. Therefore the integration over the spin-2 component of 
the metric in the present model leads to the same result as in (l67li when k\ = oj 2 + k 2 (in 
Euclidean space) is replaced by uj\ + k\ — Fi(kf) 2 — Si(/c x ) 3 . Assuming that F x and Si are 
negative constants, we obtain then 

exp {"W 

(74) 



In order to calculate the trace above, we first integrate over the frequencies 


Tr 


( S{ki + k 2 ) ) 

I * (0) 1 

r d 3 k 

V^ x 2 + k\ + 

|F 1 |(£?) 2 + |s 1 |(£?)p 

‘ 2 J 

< 2 "> 3 J¥T\ 


5(0) 


\Fi\ 


where 


Z(A) = 


d 3 k 


2 v ^ VN 

1 


( 2n ) 3 kWk 2 + A ’ 

Using dimensional regularisation, this integral becomes 

k d ~ 3 


Z(A) = 


2n d /2 ^-d r oo 

rW 2) W d Jo {k 2 + A) 1 / 2 


dk- 


which is calculated using na 


dk 


k? 


Writing d = 3 — e, we find then 


__ T(^)T(a- 

(k 2 + A) a 2 (A) Q " 1 ^r(a) 


Z < A > = if + ■ 


(75) 

(76) 

(77) 

(78) 

(79) 
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which finally leads to 


Tr 


6(ki + k 2 ) 


m 




M + k 2 + l-Fil(fcf) 2 + |^|(fcf)3 / (27T) 2 v1^ 


+ 


(80) 


where dots represent finite terms. We note here that the above result does not depend on 
|F\|, since this coupling constant controls k 4 in (1741) . which is sub-dominant, and the UV 
behaviour is dominated by the terms |S'i|/c 6 . 


Spin-0 component 

Gathering all terms which depend on h in the action (PHI) , we write the equivalent action 
in terms of the Fourier transform h of h and, after performing a Wick rotation, we obtain 


Vhexp < iS 


- 1 ? 




d 4 ki d 4 k 2 1 tn 

(27r) 4 (27r) 4 2^ " 2 


(ii) 


where 


A 


(ii) 




Vi 

P 2 (^), 


+ 2 


A 


T>! (k 


h(k x 


(81) 


- [Mi • q)M 


Vi(h 


U 1 


X 2 


8 k 2 


V 2 {k\) 

5{k\ + k 2 ) , 



(82) 


and 


kl 


KAM = Xujf + 

MM = -[kl - F 3 (kl) 2 + S 3 (kl) 3 ] 
V 2 {k x ) = ~[akl - F,(kl) 2 + S 4 (kl) 3 ] 


1 + ( 3 Fi + 8^2)^! + (3 Si + 8S 2 )(fc 1 ) 2 ^ 2 (/ci) + 2 'Diiki) 


-V 2 (ki 


(83) 


As expected, the IR limit (he. neglecting terms of the form ( k 2 ) n for n > 1) of the expression 
(1S2D leads to the expression (I70|) . and the dispersion relation for the scalar graviton is given 
by the identity (1351) . The functional integration over h gives 


exp 


2Mp L 


iP M + AMiMM 


Tr 


Mp L 
X 2 


12 M 2 L 


(pVo) + AMiMM 


Tr 


ifn) + AMl • q)ipo\ 


4 


Tr 


8{ki + k 2 ) 

K l M) 

8(k 1 + k 2 )Vi(ki) 
5(k 1 + k 2 )u\ 


, « Tr 

M 2 


(84) 


The first three traces above can be written in the following generic form 

5{ki + k 2 ) f T>i(ki) 


Tr 


K M) Vp 2 (fci) 


for n = 0,1, 2 , 


(85) 
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which, after integrating over the frequencies leads to 


Tr 


5{ki + k 2 ) ( T>i(ki) 


K\ k i) b 2 (fei) 


5(0) 


d 3 k 


2 yfX J ( 2 t r) 


m 


P2(k) 

Di$). 


2 n 


( 86 ) 


with 


k 2 


= (?V - Xu 2 ) = 


l + (3F 1 + 8F 2 )k 2 + (3S 1 + 8S 2 )(k 2 ) 2 V 2 (k) + 2 ^(fc) 


-^ 2 (fc) 


(87) 

Expanding the terms inside the square brackets in (l 86 j) . keeping only terms of at least quartic 
order in the momentum, as we did in for the tensor field, the right-hand side of eq. (l 86 j) reduces 
to 


5(0) 


a 


(n)' 


2 \ xc, 


(n) 


c, 


(n) 


with the integral Z given by eq. 04711 and, for n — 0,1, 2, 


( 4 n) /M* HL 

(4”’/v 


cf> = 

r , ( n ) — 

°6 — 

4 °‘ = —(2/i + 8/2) - 


( 88 ) 


( 89 ) 


HL ) 


2S3(/4S3 — 2 / 3 S 4 ) 


c (1) 

c 4 


— 2/4 + 


2 s 4 ( 2 si + 8S2) (/4S3 — /3S4) — 535^(2/4 + 8/2) 


C® = 


s 4 [^ s 3 s 4 ( 2 /l + 8/2) + 4s 4 (/ 4 S3 — /3S4) ( 2 <S 4 + 8S2) + 2s§(3/ 4 S3 — 2 /3S4)] 


c (0) - 


.( 1 ) _ 


2 s 2 

— (2si + 8 s 2 )-— 

s 4 


= —2s 4 — 


(2si + 8 s 2 )s 2 


J 2 ) _ 


= - 2 - 1 - 


s 4 S 4 ( 2 si + 8 S 2 ) 


Solving this integral with dimensional regularisation using the result (1481) . we obtain 


Tr 


5(ki + k 2 ) ( V 1 (k 1 ) 




5(0) 


F 


(2t)V*c“ £ 


+ • ' • , 


where dots represent finite terms. Finally we need to calculate the following trace 


(90) 


Tr 


5( ki + k 2 )uj\ 

K l {h)ki 


= m 
m 

X 


d 4 k 


LU 


( 2?r ) 4 p h\k)k 2 

r d 4 k i r 


d 4 k Q(k)/X 


( 2 ?r ) 4 % 2 J ( 2 tt ) 4 [ w 2 + g{%)/x}k 2 


(91) 
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with Q defined in f[871) . The first integral in the last line of (19T1) vanishes with dimensional 
regnlarisation. In addition, for the second integral, we integrate over the frequencies to get 


Tr 


5{k\ + k 2 )u]\ \ 5(0) f d 3 k \G(k) 


V^ik^h 


2X 3 / 2 J (2vr) 3 p 

We then expand G(k), for which only the dominant term potentially leads to a divergence 

'S(0)^Cf ] r d 3 k 


(92) 


Tr 


K\ki)ki 


2X 3 / 2 


(2vr) : 


\k\ +••• 


(93) 


with dots representing finite terms. The integral vanishes with dimensional regnlarisation, 
such that the result is finite. 


Acknowledgements The work of J. L. is supported by the National Council for Scientific 
and Technological Development (CNPq - Brazil). 


References 

[1] J. F. Donoghue, AIP Conf. Proc. 1483 (2012) 73 |arXiv:1209.3511 [gr-qc]]. 

[2] P. Horava, Phys. Rev. D 79 (2009) 084008 [arXiv:0901.3775 [hep-tli]]. 

[3] C. Charmousis, G. Niz, A. Padilla and P. M. Saffin, JHEP 0908 (2009) 070 
[arXiv:0905.2579 [hep-tli]]; D. Bias, O. Pujolas and S. Sibiryakov, Phys. Rev. Lett. 104, 
181302 (2010) [arXiv:0909.3525 [hep-th]]. 

[4] N. Arkani-Hamed, H. C. Cheng, M. Luty and J. Thaler, JHEP 0507 (2005) 029 
| hep-ph /0407034]. 

[5] M. Pospelov and Y. Shang, Phys. Rev. D 85 (2012) 105001 [arXiv: 1010.5249 [hep-tli]]. 

[6] I. Kimpton and A. Padilla, JHEP 1304, 133 (2013) [arXiv: 1301.6950 [hep-tli]]. 

[7] R. Iengo, J. G. Russo and M. Serone, JHEP 0911 (2009) 020 [arXiv:0906.3477 [hep-th]]. 

[8] J. Alexandre, J. Brister and N. Houston, Phys. Rev. D 86 (2012) 025030 
[arXiv: 1204.2246 [hep-ph]]. 

[9] J. Alexandre and J. Brister, Phys. Rev. D 88 (2013) 6, 065020 [arXiv:1307.7613 [hep-th]]. 

[10] A. Contillo, S. Rechenberger and F. Saueressig, JHEP 1312 (2013) 017 [arXiv:1309.7273 
[hep-th]]; G. D’Odorico, F. Saueressig and M. Schutten, Phys. Rev. Lett. 113 (2014) 17, 
171101 [arXiv: 1406.4366 [gr-qc]]. 

[11] G. D’ Odorico, J. W. Goossens and F. Saueressig, JHEP 1510 (2015) 126 
[arXiv: 1508.00590 [hep-th]]. 


22 










[12] M. Henneaux, C. Teitelboim and J. Zanelli, Nucl. Phys. B 332 (1990) 169. 

[13] G. W. Gibbons, S. W. Hawking and M. J. Perry, Nucl. Phys. B 138 (1978) 141; 

[14] P. O. Mazur and E. Mottola, Nucl. Phys. B 341 (1990) 187. 

[15] G. ’t Hooft, Subnucl. Ser. 40, 249 (2003). 

[16] J. Alexandre, J. R. Ellis and N. E. Mavromatos, JHEP 0612 (2006) 071 
|hep-th/06100721; J. Alexandre and N. E. Mavromatos, hep-th/0703171, 

[17] G. Leibbrandt, Rev. Mod. Phys. 47 (1975) 849. 

[18] S. Weinzierl, Mod. Phys. Lett. A 29 (2014) 15, 1430015 [arXiv: 1402.4407 [hep-ph]]. 

[19] L. H. Ryder, “Quantum Field Theory,” ISBN-9780521237642. 

[20] V. A. Kostelecky and N. Russell, Rev. Mod. Phys. 83 (2011) 11 [arXiv:0801.0287 [hep- 
ph]]. 

[21] J. Alexandre and J. Brister, Phys. Rev. D 92 (2015) 2, 024025 [arXiv:1505.01392 [hep- 
th]]. 

[22] P. Horava and C. M. Melby-Thompson, Phys. Rev. D 82 (2010) 064027 [arXiv:1007.2410 
[hep-tli]]. 

[23] G. Degrassi, S. Di Vita, J. Elias-Miro, J. R. Espinosa, G. F. Giudice, G. Isidori and 
A. Strumia, JHEP 1208, 098 (2012) [arXiv:1205.6497 [hep-ph]]. 

[24] M. Herranen, T. Markkanen, S. Nurmi and A. Rajantie, Phys. Rev. Lett. 113 (2014) 
21, 211102 [arXiv:1407.3141 [hep-ph]]. 

[25] J. Alexandre, J. Leite and N. E. Mavromatos, Phys. Rev. D 87 (2013) 12, 125029 
[arXiv:1304.7706 [hep-ph]]; J. Alexandre, J. Leite and N. E. Mavromatos, Phys. Rev. D 
90 (2014) 4, 045026 [arXiv: 1404.7429 [hep-th]]. 


23 



